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Abstract. The paper is concerned with the properties of the distance function from a closed subset 
of a liiemannian manifold, with particular attention to the set of singularities. 
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1. Introduction 

We are concerned with the properties of the singular set of the distance function from a closed 
subset of a n-dimensional, smooth and connected Riemannian manifold (M, g), with particular 
attention to its rectifiability. 

Definition 1.1. We say that a subset S of {M, g) is C^-rectifiable, with r > 1, if it can be covered 
by a countable family of embedded submanifolds of dimension (ri — 1), with the exception of 
a set of 7Y"^^ zero measure, where 7Y"^^ is the {n — l)-dimensional Hausdorff measure on M. 
We will simply say that a set is rectifiable when it is at least C^-rectifiable. 

See [18, 29] for a complete discussion of the notion of rectifiability. 

The distance function from a closed, not empty subset K of (M, g) is defined in the usual way, 

dji{x) = inf d{x,y) 
yeK 

where d is the distance on M induced by the metric tensor g. 

The singular set Sing of dx : M ^ R is the set where this function fails to be differentiable. 
Our study of the rectifiability of Sing relies on the theory of viscosity solutions of Hamilton-Jacobi 
equations. Indeed, we show that the distance function dx is a viscosity solution of the following 
problem 

f|Vu| = l inM\i^, 
I w = on dK 

and we use the property of semiconcavity shared by such solutions to obtain a rectifiability result 
for SiJig. 

Then, we investigate under which hypotheses also the closure of Sing is rectifiable. This prob- 
lem is strictly connected to the analysis of the geodesic flow on (A/, g) originating from K, hence. 
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we adapt ideas coming from the study of the cut locus of a point in a Riemannian manifold, which 
is actually the very special case when is a single point of M. 

Our results lead to the conclusion that, under some conditions on the regularity of the set K, 
the Hausdorff dimension of the closure of the singular set is at most (n — 1) and that the gradient 
of the distance function from K is locally a vector field with special bounded variation (see [3, 4, 5]). 
Moreover, we also study when the singular set shares an higher regularity and we analyse in 
detail its topological structure if M is a two-dimensional analytic surface and K an analytic 
subset. 

The study of the distance fimction and of the associate eikonal equation | Vu| = 1 is a special 
example of a connection which can be extended to a large class of stationary Hamilton-Jacobi 
equations. In the last section we discuss some problems about the structure of the singular set 
of more general viscosity solutions, suggested by some geometric results for the cut locus of a 
point. 

2. Stationary Hamilton-Jacobi Equations on Manifolds 

Let M be a smooth and connected, n-dimensional differentiable manifold. 
We consider the following Hamilton-Jacobi problem in c M, 

{ll{x,du{x),u{x)) = inU, 
u = uq on dfl 

where H : T*il x R ^ M and T* denotes the cotangent bundle. 

Definition 2.1. Given a continuous fimction u : —^R and a point x G M, the superdijferential of 
M at a; is the subset of T*M defined by 

d'^u{x) = ^dip(x) I if G C^{M),(p(x) — u{x) = mins/? — w| • 

Similarly, the set 

d~u{x) = < dil){x) I il) £ C^iM),il){x) — u(x) = maxip — u} 
is called the subdifferential of u at y. 

Notice that it is equivalent to replace the max (min) on all M with the maximum (minimum) in 
an open neighborhood of x in M. 

It is easy to see that d'^u{x) and d~u{x) are both nonempty if and orvly if u is differentiable at 
X G M.ln this case we have 

d~^u{x) = d~u{x) = {du{x)} . 
We list here without proof some of the standard properties of the sub and superdifferentials 

which will be needed later. 

Proposition 2.2. If ip : N ^ M is a map between the smooth manifolds N and M which is around 
X G N, then 

d~^{u o tp){x) D d~^u{'tjj{x)) o dil){x) = {v o dil){x) \ v G d'^u{ip{x))} . 

Iftp is a local diffeomorphism near x, the inclusion becomes an equality. An analogous statement holds for 
d-. 

Proposition 2.3. If9:lS.—flS.isa function such that 0{u{x)) > 0, then 

d+{9 o u){x) D d9{u{x)) o d+u{x) = {de{u{x)) ov\vG d+u{x)} , 

similarly for . If9{u{x)) > then the inclusion is an equality. 

For a locally Lipschitz fimction m on a Riemannian manifold {M,g), d^u{x) and d~u{x) are 
compact convex sets, almost ever5rwhere coinciding with the differential of the function u, by 
Rademacher's Theorem. 

For a generic continuous function u we prove in the next proposition that d^u{x) and d~u{x) are 
not empty in a dense subset. 
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Proposition 2.4. Let u : ft ^Rbea continuous function on an open subset O ofM. Then the subdijfer- 
ential d^u{x) (the superdifferential d^u{x)) is not empty for every x in a dense subset ofQ.. 

Proof. It is always possible to endow M with a Riemannian structure giving a metric d{- , •) on M 
which generates the same topology. 

Consider a generic point y E and a geodesic ball B contained in with center y. If the ball B 
is small enough, the function x i-^ d^{x, y) is smooth in B. Taking a large positive constant A, 
the function Fa (a;) = u{x) + Ad'^{x, y) has a local minimum at a point xa in the interior of B. At 
XA the subdifferential of the function Fa must contain the origin of T*^M, hence, being y) 
differentiable in the ball B, the differential of —d'^{x, y) at xa belongs to d~u{xA)- As the point 
y and the ball B were arbitrarily chosen, the set of points where the subdifferential of u is not 
empty is dense in fi. 

The same argument holds for the superdifferential of u, considering the function —u. □ 

Now we introduce the notion of semiconcavity which will play a central role in the first part 
of the paper. 

Definition 2.5. Given an open set Q. C M", a continuous function w : f2 ^ M is called locally 
semiconcave if, for any open convex set Q' c ^ with compact closure in f2, there exists a constant 
C such that one of the following three equivalent conditions is satisfied, 
1. Va;, h with x,x + h, x — h e il', 



2. u{x) — C\x\'^ is a concave ftinction in Q', 

3. D'^u < 2C Id in O', as distributions {Id is the n x n identity matrix). 

In order to give a meaning to the concept of semiconcavity when the ambient space is a differ- 
entiable manifold M, we analyse the stability of this property under composition with maps. 

Proposition 2.6. Let fl and Q' two open subsets ofW\ Ifu:^^Risa Lipschitz function such that 
u{x) — C |a;p is concave and tp : ^' ^ ^ is a C"^ function with bounded first and second derivatives, then 
uoip : fl' -^Ris a Lipschitz function and u o tp^y) — C"|t/p is concave, for a suitable constant C . 

The proof is straightforward. 
Then, the following definition is well-posed. 

Definition 2.7. A continuous fimction m : M — > M is called locally semiconcave if, for any local 
chart : K" ^ f2 c M, the function uoipis locally semiconcave in M". 

The importance of semiconcave functions in connection with the generalized differentials is 

expressed by the following proposition (see [12]). 

Proposition 2.8. Let the function u : M ^ Rbe locally semiconcave, then the superdifferential d^u is 
not empty at each point, moreover, is upper semicontinuous, namely 

Xk X, Vk ^ V, Vk e d'^u{xk) V e d'^u{x) . 

In particular, if the differential du exists at every point ofQ e M, then u e C^(f2). 

Now we introduce the definition of viscosity solution. 
Let ri be an open subset of M and H, called Hamiltonian fimction, a continuous real function on 
T*f2 X R. We are interested in the following Hamilton-Jacobi problem 



Definition 2.9. We say that a continuous function u is a viscosity solution of equation (2.1) if for 
every x gQ, 



If only the first condition is satisfied (resp. the second) u is called a viscosity subsolution (resp. a 
viscosity super solution). 



u{x + /i) + u{x -h)- 2u{x) < 2C\h\^ 



(2.1) 



H(a;, du{x), u{x)) = 



in fl . 



(2.2) 
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If fl' is an open subset of another smooth differentiable manifold and ?/> : fi' ^ f2 is a 
local diffeomorphism, we define the pull-back of the Hamiltonian function ip*B. : T*fi' x R ^ R 
by 

ip*'R{y, V, r) = H(V'(?/), v o dtpiyy^r) . 
Taking into account Proposition 2.2, the following statement is obvious. 

Proposition 2.10. If u is a viscosity solution of R = in c M and tp : Q' ^ Q is a local 
diffeomorphism, then uoipisa viscosity solution o/V'*H = in Ct' c A''. 

3. The Distance Function from a Subset of a Manifold 

From now on, (M, g) will be a smooth, connected and complete Riemannian manifold without 

boundary, of dimension n. 

We will study the distance fimction dK from a closed and not empty subset K of M (for tech- 
nical reasons, sometimes we consider also its square dj^ix) : M ^ R). 

The distance between two points x and y or from the point x to the set K is defined as the infi- 
mum of the lengths of the curves starting at x and ending at y, or on K, respectively As M is 
complete, by the Theorem of Hopf-Rinow, such infimum is reached by at least one curve which 
will be a smooth geodesic. 

The distance from the set i^T is a continuous function on M but in general it is not everywhere 
differentiable, for instance, if the manifold M is compact, the distance function from any proper 
subset will be singular at the points of absolute maximum. This section deals precisely with the 
set where the gradient of dx does not exist. 

In the following we will denote the distance between two points x,y G M with d{x, y) and 
the exponential map of {M,g) with Exp : TM xR M. For simplicity, we will write \v\ for the 
modulus of a vector v S TM, defined as \Jg{v, v). 

Moreover, we will take often into account the identification between the differential du and the 
gradient Vu of a fimction, induced by the scalar product g. 

Theorem 3.1. The distance function dx is the unique viscosity solution of the following Hamilton-Jacobi 
problem 



(3.1) 



|Vm|-1=0 inM\K, 
u = on K 



in the class of continuous functions bounded from below. 
The function dj^/2 is the unique viscosity solution of 



(3.2) 



|Vu|-2u = inM, 
u = on K 



in the class of continuous functions on M such that their zero set is K. 

Remark 3.2. The restriction to lower bounded functions is necessary, and — are both vis- 
cosity solutions of Problem (3.1) with M = R" and K = {0}. Moreover, the completeness of M 
plays an important role here, if M is the open unit ball of R" the same example shows that the 
uniqueness does not hold. 

Notice also that every function where ff is a closed subset of M with if D iC, is a viscosity 
solution of Problem (3.2), equal to zero on K. 

Proof. Notice that dxix) is the minimum time t > for any curve 7 to reach a point ^{t) G K, 
subject to the conditions 7(0) = and I7I < 1; dx is then the value function of a "minimum time 
problem"; this proves that dx is also a viscosity solution of Problem (3.1), by well known results 
(see for example Proposition 2.3, Chapter IV in [13]). Then we show that the fimction d\/2 is a 
solution of Problem (3.2). 

First of all, notice that the distance function from if is a 1-Lipschitz function, hence d\ is locally 
Lipschitz. 
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As dK is 1-Lipschitz, at every point of K the function dj^- is differentiable and its differential is 
zero. Hence, the definition of viscosity solution holds for points belonging to K. In order to 
prove the theorem, it is then sufficient to test conditions (2.2) on the generalized differentials at 
the points of the open set M\K. 

Since d|-/2 is positive in M \ i^, applying Proposition 2.3 with the function 6{t) = \f2±, we see 
that the function d|-/2 is a viscosity solution of 

V\/2^ \/2^/ 
in M \ if. Being there positive, it also solves 

6i(Vm, Vm) - 2u = 

\r\M\K. This fact together with the previous remark about the behavior of at the points of 
K gives the claim. 

Suppose now that u is a viscosity solution of Problem (3.1) then, u is also a solution of 

('|Vu|-l = vsxM\K, 
[u = ox\.K. 

As in the work of Kruzhkov [23], we consider the fimction v = — e~" which, by Proposition 2.3, 
turns out to be a viscosity solution of 

j\Vv\+v = in M\K, 
yv = —1 on K 

moreover, \v\ < e~'"^". 

We establish an uniqueness result for this last problem in the class of boimded functions v, which 
clearly implies the first uniqueness result. We remark that the proof is based on similar ones 
in [14, 15, 20]. 

We argue by contradiction, suppose that u and v are two boimded solutions of (3), \u\, \v\ < C, 
and that at a point x we have u{x) > 2e + v{x) with e > 0. 
Let b{x, y) : M x M ^ M be a smooth function satisfying 

• 6 > 

. \S/^h{x,y)\,\\/yh{x,y)\ <2 

• supmxm \d{x,y) - h{x,y)\ < oo 

such a function can be obtained smoothing the distance function in M x M. 

We fix a point xo in K and we define the smooth function B{x) = xq)"^. By the properties of 

b and the boundedness of u and v, the following function : M x Af ^ M 

y) = u{x) - v{y) - Xd{x, yf - 5 B{x) - 5 B{y) 

has a maximum at a point x, y (dependent on the positive parameters 6 and A) and such maxi- 
mum ^{x, y) is less than 2C. Hence, the fimction 

(3.3) x ^ [v{y) + Xd{x, yf + 5 B{x) + 5B{y)] - u{x) 
has a minimirai at x while 

(3.4) y ^ [u{x) - \d{x, yf - 5B{x) - 5B{y)] - v{y) 

has a maximum at y. 
If 25 < e/B{x) then 

^{x,y) > *(S,S) > 2e - 25B{x) > e 

hence, we get 

(3.5) 5B{x) + 8B{y) + Xd{x, yf+e< u{x) - v{y) < 2C . 

This shows that, for a fixed 5, the pair x, y is contained in a bounded set and, if A goes to +oc the 
distance between x and y goes to zero. Possibly passing to a subsequence for A going to infinity, 
X and y converge to a common limit point z which cannot belong to K, otherwise we would get 
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e < u{z) - v{z) = 0, thus, for some A large enough also x and y do not belong to K. 
As the function tf {x, y) is smooth in x c M x M, where is a small geodesic ball around 
z, choosing a suitable A large enough we can differentiate the functions inside the square brackets 
in equations (3.3) and (3.4) obtaining 

v = 5 VB{x) + \V:^(f{x, y) e d+u{x) , 

w = -SVB{y) - \Vy(f{x,y) G d~v{y) . 

By Definition 2.9 we have that + u{x) < and \w\ + v{y) > 0, hence 

u{x) — v{y) + \v\ — |w| < . 

Moreover, 

\v\ - \w\ = \SVB{v) + XV,d''{x,y)\ - \dVB{y)+\Vyd^ix,y)\ 

> \XV,d\x,y)\ - \XWyd^{x,y)\ - \5VB{y)\ - \5VB{x)\ 

= 2Xd{x,y) \V,d{x,y)\ - 2Xd{x,y) \Vyd{x,y)\ - \SVB{y)\ - \dVB{x)\ 

= 2Xd{x,y) - 2Xd{x,y) - \dVB{y)\ - \SVB{x)\ 

= -\SVB{y)\-\SVB{x)\ 

which implies, 

u{x) - v{y) - d\VB{y)\ - 6\VB{x)\ < 0. 

Finally, we have that 

6\VB{x)\ = 26\b{x,xo)yb{x,xo)\ < 4,6y/B{x) 
and using the estimate 5B{x) < 2C which follows from equation (3.5), 

5\VB{x)\ < d,V25C < e/4 
if 6 was chosen small enough. Holding the same for y, we conclude that 

u{x) - v{y) - £/2 < 

which is in contradiction with the fact that ti,{x) — v{y) > e. 

About the second xmiqueness claim, if u is a continuous viscosity solution of Problem (3.2) 
then, by Proposition 2.4 the superdifferential of u is not empty in a dense subset of M \ if, hence, 
directly by the equation and by continuity, u is non negative. By the h5qDothesis on its zero set 
we conclude that u is positive mall M\K. Composing u with the function 1 1—^ \/2i, we see that 
\/2u is a positive, continuous viscosity solution of Problem (3.1), then it must coincide with dK, 
by the previous result. It follows that u = d%/2. □ 

We now study the rectifiability of the singular set of dK, 

Sing = {a; e M I is not differentiable at a;} . 

Remark 3.3. In this definition we used the squared distance function instead of the distance in 
order to avoid to consider also the points of the boimdary of K, which are singular for dK- It is 
easy to see that outside K the distance and its square have the same regularity. 

Proposition 3.4. The function dK is locally semiconcave in M\K. 

Proof. The distance function dK is a viscosity solution of H = in M \ if , where the Hamiltonian 
function is given by H(a:, u, t) = |up — 1. We choose a smooth local chart -0 : ^ C M and we 
define v = dK ° i^, which is a locally Lipschitz function and, by Proposition 2.10, it is a viscosity 
solution of V'*H = 0. 

The pull-back of the Hamiltorvian function on R" takes the form 

V'*H(?/,w,s) = s'^(y)(rfV(w),#(w)) - 1 = gij{y)wiWj - 1 

for (y, ui, s) e M" X R" X R and where {y) are the components of the metric tensor of M in local 
coordinates. 

Since the matrix gij{y) is positive definite ip*'H.{y,w,s) is locally uniformly convex in w, hence. 



HAMILTON-JACOBI EQUATIONS AND DISTANCE FUNCTIONS 



7 



by Theorem 5.3 of [24], it follows that v = dx is locally semiconcave in M". Recalling Defini- 
tion 2.7, this means that cIk is locally semiconcave in M\K. □ 

The semiconcavity of cIk allows us to work with the superdifferentials when the gradient does 
not exist. Indeed, notice that the points of Sing are precisely those where the superdifferential is 
not a singleton, then, the following result is a straightforward consequence of Proposition 2.8. 

Proposition 3.5. The function d\ is of class in the open set M \ Sing and dx is in M \ 
(KuSing). 

Remark 3.6. The semiconcavity property also gives information about the propagation of singu- 
larities and the relations between the structure of the superdifferential at a point x and the set of 
minimal geodesies from x to K (see [2, 6]). 

The set Ext{d'^d\{x)/2) of extremal points of the (convex) superdifferential set of c?|./2 at x is 
in one-to-one correspondence with the family Q {x) of minimal geodesies from x to K. Precisely 
Q{x) is described by 

(3.6) g{x) = {Exp(a;, -v, ■) : [0, 1] ^ M | for v e Ext{d+d]i{x)/2)} . 

The set of points of K at minimum distance from x are given by Exp(x, —v,l) for v in the set of 
extremal points of the superdifferential set of d\ /2 at x. As a particular case we have that if the 
function d\ is differentiable at x, then the point of K closest to x is uniquely determined and 
given by Exp(x, -V(i|(x)/2, 1). 

Finally, notice that Sing is the set of points x such that the distance d{x, K) is realized by more 
than one minimal geodesic between x and K. 

The rectifiability of Sing now follows. 

Proposition 3.7. The set Sing is C^-rectifiable. 

Proof. By a result proved in [1], the singular set of a locally semiconcave function in an open set 
of R" is C^-rectifiable. We take a countable family of local charts ^'i : and consider the 

functions dx oipi. These functions are locally semiconcave in M" with singular sets Singi, hence, 
by the relation 

oo 

Sing C [J tpiiSingi) 

i=l 

we get the thesis. □ 

The same statement does not hold for the closure of Sing, for a generic closed set. We describe 
now a counterexample showing that indeed the set Sing is not rectifiable for a set K of class C^'^ 
only. 

We look for a convex open set O with a C^' ^ boundary in such that the closure of the singular 
set Sing of the distance function from its boundary has nonzero Lebesgue measure, hence it is 
not rectifiable. 

We start with a Cantor-like set C C S^, closed with empty interior in S\ with no isolated points 
and positive Hausdorff measure. Such a set can be constructed as follows 

oo 

C = S'\[jli 

i=l 

where {/»} is a coimtable family of open disjoint connected arcs on whose middle points are 
Pi e and such that the sum of their lengths is less than 27r. 

We claim that every point of C is a limit point of the sequence {pi}- If p E C there must be a 
sequence of arcs li- arbitrarily close to p, since the arcs are countable and the sum of their lengths 
is boimded by 27r we have that they shrink when j goes to infinity, hence pi- — > p. 
We define an open convex set Q' as the intersection of the open halfplanes, containing the origin 
of M^, determined by the tangent lines to at the points of C, see the following figure. 
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Let us take an arc I with middle point p, bounded by P, Q e C and consider the associate quadri- 
lateral OPVQ. If the point x is inside the open triangle OPV it is clear that the point of dVt' 
closest to X belongs to the segment PV and it is unique (a? in the figure). Hence, for such points 
the distance from the boundary of Q.' coincide with the distance from the segment PV. 
Applying the same argument to the open triangle OQV , we see that the segment OV consists of 
singular points of dga', moreover, the segment OV intersects §^ at the point p. 
It follows that the imion S of the segments from the middle points Pi to the origin coincides with 
Sing for dgw {x). Being pi dense in C, the closure of S contains AC C for every A e [0, 1]. As C 
has Ti} positive measure, the Lebesgue measure of S is positive, by Fubini's Theorem. 
Now let Q. be the set of points of whose distance from the convex Vt' is less than 1. It is 
immediate to check that 

ddn {x) = ddo,' {x) + 1 Mx eQ.' 
hence for every x in the unit ball. 

So the closure of Sing for the distance function from the botmdary of O (or from the comple- 
mentary set of fi in M?) has positive Lebesgue measure, moreover, by the properties of convex 
bodies, the boimdary of VI is of class at least C^'^. Since the Lebesgue measure of Sing is positive 
it cannot be rectifiable. 

Remark 3.8. In the next section we will show that if the boimdary of K is of class at least 
then also the closure of Sing is rectifiable. To our knowledge it is unknown even in if the 
gap between such result and the previous counterexample can be filled, that is, if the (or 
maybe C^'^) regularity of the boimdary of K is enough to get the rectifiability of the closure of 
the singular set. 

4. Rectifiability of the Closure of the Singular Set 

In this section we are going to show that an higher regularity of the set K implies the recti- 
fiability also of the closure of the singular set. Moreover, we determine a relation between the 
regularity of K and of the hypersurfaces covering Sing. 

In all this section K is a fc-dimensional embedded C" submanifold of M without boundary, 

with < A: < n — 1 (the case fc = n is trivial) and r > 2. 

Let UK be the unit normal bundle of in M, we denote with F : UK x [0, +oo) M the 
restriction of the exponential map of M to UK x [0, -l-oo). Since K is C, UK is a manifold of 



HAMILTON-JACOBI EQUATIONS AND DISTANCE FUNCTIONS 



9 



class C" ^ and being Exp a smooth map, F{q, v, t) and all its derivatives in the t variable are C" ^ 
functions. 

Remark 4.1. The case when K is the closure of an open set of M with smooth boundary can be 
reduced to our case. Indeed, if a; e M \ K, then the minimal geodesic from x io K ends in 
dK without touching the interior points, hence to study dx we can simply consider the distance 
function from dK in every open cormected component of M \ K. However, in this case, some 
results like the higher smoothness of d\ at the points of dK expressed by Proposition 4.2, could 
be lost. 

The behavior of near K is weU known. 

Proposition 4.2. Thereexiste > OandanopenneighborhoodCtofK in M suchthatthemap F\ukx{o,£) '■ 

UK X {0,e) —^ il\K is a C""~^ dijfeomorphism. 

Moreover, 

• for every point in O there is an unique point of minimum distance in K (the unique projection 
property holds for K in Q.), 

• the distance function dx is C" in fl\K, 

• the squared distance function d% is in Q. 

Remark 4.3. It can be proved that K has to be at least (7^'^ in order to share the unique projection 
property in a neighborhood, in such case the squared distance function also turns out to be of 

class Ci'i. 

See [16, 17] for a detailed discussion of the relation between the regularity of K and of dK- 

In order to study what happens /ar from K we have to analyse the sets of points where the 

unique projection property fails or F is not a local diffeomorphism. From a topological point of 
view, the problem is naturally cormected with the study of the singularities of maps between 
Euclidean spaces. For instance, when K is a single point of M the singular sets were shown to be 
related to the classes of singularities considered by the Theory of Catastrophes, see [11]. 

Consider the geodesic curve 1 1-^ F{q, v, t) for t G [0, to] {{q, v) G UK is fixed), for small values 
of to it is the imique minimizer of the length functional between its end point p = F{q,v,to) and 
K but for large to > it could cease to be minimal. Hence, there exists a value a (possibly +oc) 
such that this geodesic is minimal between q and F{q, v, t) for every t < a, but not on any larger 
interval. If a{q, v) < +oo, we say that the point F{q, v, a{q, v)) is the cut point of the geodesic 
F{q, V, t) and we define the following set, 

Vk = {{q,v,t) eUKxM.+ \t< (j{q,v)] . 

Notice that the set FIVk) clearly contains Vl\K, where Q is the open neighborhood of Proposi- 
tion 4.2. 

Definition 4.4. The set of points F{q, v, a{q, v)) for (g, v) G UK with a{q, v) < +oo is called the 
cut locus of K, we denote it with Cut{K). 

The reasons why a geodesic ceases to be minimal are explained by the following proposition 

(see [19, 28]). 

Proposition 4.5. If for a geodesic F{q, v, t) we have a{q, v) < +oo, at least one of the following two non 
exclusive conditions is satisfied: 

1. at the point p ~ F{q, v, a{q, v)) there arrives another minimal geodesic from K, 

1. the differential dF{q. v, a{q, u)) is not invertible. 
Conversely, if at least one of these conditions is satisfied the geodesic F{q, v, t) cannot be minimal on an 
interval larger that [0, a{q, v)]. 

Notice that, by Remark 3.6, if condition 1 above is satisfied then the point p belongs to Sing, 
while Sing is clearly included in Cut{K). 

Then, we consider the following two subsets of Cut{K), 
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• Sing is the set of points p = F{q, v, a{q, v)) where more than one minimal geodesic from K 
arrives, 

• J is the set of points p = F{q,v,cr{q,v)) such that the differential dF{q,v,a{q,v)) is not 
invertible. 

We call J the locus of optimal focal points. 
Clearly Sing U J = Cut{K). 

The next proposition (see [19, 28]) establishes the connection between Cut{K) and the distance 
function from K. 

Proposition 4.6. The following statements hold, 

1. Cut{K) = Sing, that is, the cut locus of K is closed in M and Sing is a dense subset. 

2. The set Vk is open in UK x R+. 

3. The map a : UK R+ is continuous. 

4. The map F is a C"""^ bijection between Vk and M\{K U Cut{K)) with a C"""^ inverse. 

5. The cut locus Cut{K) is equal to F^OVk) where the boundary is considered in the ambient space 

UK X M+. 

6. The set M \ Cut{K) can be continuously retracted on K, and, ifa{q, v) < +oo for every {q, v) G 
UK, the set M\K can be continuously retracted on the cut locus Cut{K). 

7. The open set M \ Cut (K) has the unique projection property, moreover the squared distance function 
d\ is of class C in it. The distance dx is C in M \ {Cut{K) U K). 

By the point 1, the closure of Sing is precisely the cut locus of K, which is the tmion of Sing 
and J. We study then the rectifiability of these two sets separately, as the rectifiability property is 
clearly stable under countable union (see Definition 1.1). 

We have seen in Proposition 3.7 that the set Sing is always C^-rectifiable. We partially improve 
this result when K is more regular. 

Proposition 4.7. IfK is of class C" with r >3the set Sing \J isa C^-rectifiable subset ofM. 

We need a preliminary lemma. 

Lemma 4.8. If there are infinitely many minimal geodesies from p to K, then p is an optimal focal point. 

Proof. If F{qi, Vi, a{qi,Vi)) = p for infinite distinct geodesies F{qi, Vi, t), then a{qi,Vi) = d{qi,p) = 
dxip) hence, by compactness, we may assume that {qi,Vi) — > {q,v) for some {q,v) e UK. It 
follows (by the semicontinuity of the length functional) that F{q,v, t) is a minimal geodesic for p 
and that dF{q, v, t) is singular since F is not locally injective near {q, v, t). □ 

Proof of Proposition 4.7. Let p he a point in Sing \ J. We know that the number of minimal 
geodesies F{qi, Vi, t) from p = F(qi,Vi, a(qi,Vi)) to K is finite by the lemma above and greater 
than one, by the singularity at p. Moreover, the differential dF{qi,Vi, <7{qi, Vi)) is invertible for ev- 
ery i, then F is locally invertible in the neighborhood of every point (gj, Vi, a{qi, Vi)) G UK x IR+. 
Let Ui be disjoint open neighborhoods of {qi, Vi, a{qi, Vi)) such that F restricted to every C/, is a 
C""^^ diffeomorphism with its image. We can also suppose that F{Ut) = U where U is an open 
neighborhood of p in M. We define the functions di : U ^ M+ given by 

d,(a;) =7rR+ {F-\x)nU,) 

where 7rR+ : UK x R+ — » IR+ denotes the projection on the second factor. Applying Propo- 
sition 4.2 to a small neighborhood of every qi, we can see that all the fimctions di are of class 
C. 

The singular set SingHU is clearly contained in the imion of the sets Sij = {r G U \ di{x) = dj (x)} 
for i j. We now prove that such sets are locally C" hj^ersurfaces. By the implicit function the- 
orem, it is sufficient to show that Vdi{x) - Vdj (x) ^ at the points of Sij. If r = F{si,Wi,t) = 
F{sj,Wj,t), for (sj, Wi, t) G Ui and {sj,Wj,t) G Uj then 

dF dF 
Vdi (x) =Vdj{x) —{si,Wi,t) = — {sj , Wj , t) . 
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Hence, by the imiqueness of the geodesic from r with a certain initial velocity vector, we would 
get {si, Wi) = {sj,Wj), contradicting the h5^othesis that Ui and Uj are disjoint. □ 

The next proposition gives the rectifiability of the set J of optimal focal points which implies 
the rectifiability of the cut locus of K, by the previous discussion. 

Proposition 4.9. IfK is of class C with r > 3, then the set J is C^'^^-rectifiable. 

Theorem 4.10. IfK is of class with r > 3, the closure of Sing, that is, the cut locus of K, is C""~^- 
rectifiahle. 

Proof of Proposition 4.9. We introduce the set J of the first focal points as follows. Let F{q, v, t) be 
a geodesic from K with {q, v) G UK and t G M+, considering the first value t = c{q, v) such that 
dF{q, V, t) is not invertible or setting c{q, v) = +oo if dF(q, v, t) is invertible for every t e M+, we 
define the map c : UK ^ M+ U {+oo}. If c{q, v) < +oo we say that F{q, v, c{q, v)) is the first focal 
point along the geodesic F{q, v, t). 

We consider the following set of points GinUK x M+ where dF is not invertible 

G = {{q,v,c{q,v)) e UK x m.+ \c{q,v) < +oo] 
and we call J = F{G) locus of the first focal points of K. 

By Proposition 4.5, we have that J D J, the set of optimal focal points, hence, it sufficient to show 
that the set J is rectifiable to prove the same for ./ and conclude the proof. 

At the points of the set G the rank of dF is at most (n — 1). We split G in two subsets, 

Gi = {{q, V, c{q, v)) e G \ Rank dF{q, v, c{q, v)) — n — 1} 
G2 = {{q, V, c{q, v)) € G \ Rank dF{q, v, c{q, v)) <n-\} . 

The following version of Sard's Theorem can be found in the book of Federer [18, Theo- 
rem 3.4.3]. 

Lemma 4.11. Let F : M" ^ M" be a map of class for some / > 1. 
If we set for any fc e {0, 1, . . . , n — 1} 

Ah = {xeW\ Rank dF{x) < k] , 

then the Hausdorff measure ofF{Ak) is zero. 

Considering local charts of UK x R+ and M and applying this lemma to our map F which is 
of class C^-i we get that 7Y"-2+2/(r-i) {F{G2)) = 0, which implies H"-^ {F{G2)) = as r > 3. 

Now we show that the set Gi is locally a C"""^ hj^ersurface in J7i^ x IR+ using the implicit 
function theorem, this proves that F{Gi) is a C"'~^-rectifiable set in M and so J. 

Let {q, V, c{q, w)) be a point in Gi, by the lower semicontinuity of the rank, choosing a small 
neighborhood B of (g, v, c{q, v))in.UK x IR+ we can suppose that there are no points of G2 in B. 
The points of Gi in S can be characterized as the zero set of the determinant det dF{q, v, t) which 
is of class C"^. 
We claim that 

..^ d det dF , , , , , „ 

(4.1) ^ {q,v,c{q,v))i^Q 

at the points of G\ n B. By the implicit function theorem, this fact implies that Gi fl B is a regular 

(n — 1) -dimensional submanifold of class C"^ and we are done. 

Let V be the covariant derivative of (M, g). Any vector in T(^g^y)UK can be represented by 
the velocity vector {w,u) = Vs(g(s), v(s))|s=o at s = of a G^~^ curve {q{s),v{s)) in UK, with 
(g(0), v{0)) = {q, v). Such a curve is given by a C curve q{s) in K with a C^~^ unit vector field 
v{s) defined along q{s) and normal to K. It is then clear that w = q'{0) belongs to the tangent 
space to K at q. 
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We set u{s) = Vsu(.s) and m(0) = u. Suppose that z{s) is an arbitrary vector field along q{s) 
tangent to K with z(0) = z, then, by the orthogonality of v{s) and z{s), we have 

= v{s))] = g{z{s), u{s)) + g{Vszis) , v{s)) 

and at the point s = we get 

= g{z, u) + giVszis), v)\s=o ■ 

Introducing the shape operator Ay : TgK —^ TgK of K at q, relative to the urvit vector v, we can 
rewrite this equation as 

g{z, u) +g{A^z, w) =0 
and, by the symmetry property of the shape operator, 

g{z, u) +g{AyW, z) = . 

Hence, as z can be chosen arbitrarily in TqK, we obtain that u + AyW € NqK. Notice also that, 
since u is a unit vector, differentiating g{v{s), v{s)) ~ 1 we obtain g{u, v) — 0. 
Resuming, the tangent space to UK at the point (g, v) is represented by the pairs of vectors 
{w, u) e TgM X TgM such that 

(4.2) w e TgK, u + AyW e NgK and g{u, v) = Q. 

Consider now a vector {w, u) e U K, and the vector field 

X{t) = duKF{q,v,t){w,u) 
along the normal geodesic 7(f) = F{q, v, t) with unit velocity vector 

i{t)=dtF{q,v,t), 

where duK denotes the partial derivative with respect to the variable (g, v) € UK. 

The field X is a Jacobi field along the geodesic 7(t), that is, it satisfies the following relations 

X{0) = w, X'{0) ^ u, 

^^■^^ X"{t)+R{X{t),i{t)W{t)=Q, 

where R is the Riemann curvature operator of M and we adopted the convention of denoting 

with T' the covariant derivative along the geodesic 7(t) of any vector or tensor field T. 

We take a basis {(wj, Uj)}, for i = 1, . . . , n — 1, of the tangent space Tf^g^^^UK and we construct an 

n-vector u) along 7 as follows, 

=Xi(t) A---AX„(t) 

where the fields Xi{t) = duKF{q, v, t){wi,Ui) are Jacobi fields and Xn{t) = 7'(t). Notice that the 
relation 

X:'Jt)+R{X^{t),-f'{t)W{t) = 

is satisfied by the field X„, since 7" — 0. As {{wi,Ui)} is a basis of Tf^g^^^UK, proving equa- 
tion (4.1) is equivalent to show that w'{c{q, v)) ^0. 

We argue by contradiction, by the Gauss Lemma (see [19, Chapter 3, Section E]) we have that 
Xn{t) is orthogonal to Xi{t) for every i = 1, ... ,n — 1, so we can suppose that the nonzero 
vector {wi,ui) is the generator of the kernel of duKF{q,v,c{q,v)). Hence, Xi{c{q,v)) = and 
X2(c{q, ?;)), . . . , Xn{c{q, v)) generate a subspace of dimension (n — 1), by the assumption on the 
rank of dF{q, v, c{q, v)). 
Computing w'{c{q, v)) we get 

uj'{c{q, v)) = X[{c{q, v)) A X2{c{q, v)) A ■ ■ ■ A Xr^{c{q, v)) 

by linearity and the h5^othesis that Xi{c{q, v)) =0. 

To conclude we need only to show that X[{c{q,v)) carmot belong to the (n — l)-dimensional 
subspace generated by X2{c{q, v)),. . . , Xn{c{q, v)). 

First we show that X[{t) is orthogonal to Xn{t) for every t. The derivative of the function h{t) = 
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g{X[{t), ^'{t)) is zero by 7" — and the last equation in (4.3), moreover the last condition of 
equation (4.2) shows that /i(0) = 0, hence the function h is identically zero. 
Consider now the function f{t) given by 

f{t)=g{X[{t),Xm-9{X,{t),X[{t)). 

We have /(O) = g{ui , Wi)— g{wi , Ui) so, using the second relation in (4.2) and taking into account 
that Wi e TqK it i < n — 1, we obtain /(O) = g{wi, AyWi) — g{AyWi, Wi) which is zero since the 
shape operator is synvmetric. Moreover, 

fit) =ff(X('(t), X,{t))-g{X,{t), Xl'it)) 

= R{Xi{t),'y'{t),j'{t),X,{t)) - R{Xr{t),i{t),i{t),Xi{t)) = 

by the properties of the curvature tensor. 

Hence, the function / is identically zero and f{c{q, v)) = gives 

giX'Mq,v)), Xi{ciq,v)))=i) for 7. = 2, . . . , n - 1 , 

so X[{c{q, v)) is orthogonal to every vector Xi{c{q, v)) and cannot belong to the subspace spanned 
by {Xi{c{q, v))} for i — 2, . . . ,n. li X'{c{q, v)) would be zero, then by the differential relation (4.3), 
we would get X(t) — X'(t) = for every t and in particular for t = 0, that is, {wi, ui) = (0, 0) 
contradicting the initial hypothesis. □ 

By standard arguments of geometric measure theory (see [18]), the rectifiability of the cut locus 
has the following immediate consequence. 

Corollary 4.12. The Hausdorff dimension of Sing {the cut locus of K) is at most {n — 1). 

To explain another consequence we need to introduce briefly the theory of functions with 
bounded variation, see [18, 29] for details. We say that a function u : R" R™ is a fimction with 
locally boimded variation u G BMoc/ if its distributional derivative Du is a Radon measure. Such 
notion can be easily extended to maps between manifolds using smooth local charts. 
A standard result says that the derivative of a locally semiconcave function stays in BVioc, in 
view of Proposition 3.4 this implies that the vector field Vdx belongs to SVJoc in the open set 
M\K. 

Now we define the subspace of BVioc of fimctions (or vector fields, as before) with locally 
special bounded variation SBVioc (see [3, 4, 5]). 

The Radon measure representing the distributional derivative Du of a function w : R" — !■ with 
locally bounded variation can be always uniquely separated in three mutually singular measures 

Du = D^u + Ju + Cu 

where the first term is the part absolutely continuous with respect to the Lebesgue measure 
Ju is a measure concentrated on a (n — l)-rectifiable set and Cu (called the Cantor part) is a 
measure which does not charge the subsets of Hausdorff dimension (n — 1). 
The space SBVioc is defined as the class of fimctions u e BVioc such that Cu = 0, that is, the 
Cantor part of the distributional derivative of u is zero. 

Corollary 4.13. IfK is of class C with r > 3, the vector field Vdx belongs to the space SBV\oc{M \ K) 
of vector fields with locally special bounded variation. 

Proof Being the cut locus rectifiable, hence of Hausdorff dimension (n — 1), the Cantor part of 

the distributional derivative of Vdx cannot be concentrated on it, so it must be concentrated in 
the open set M\{KU Cut{K)). By point 7 of Proposition 4.6, the field VUk belongs to in 
M \{K U Cut{K)) then, by the hypotheses, it is at least C^, hence its distributional derivative 
coincides with the product of the classical derivative with the Lebesgue measure, this shows that 
Cu(M \ (K U Cut{K)) = 0. These two facts together prove that Vdn belongs to SBVioc{M \ 
K). □ 
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A very important particular case is of our discussion is K ~ {p}. The cut locus of a point p 
in M arises naturally in various geometric problems, its definition is due to Poincare [27] and 
its properties were studied by many authors, see for instance [10, 11, 21, 22, 25, 26, 30, 31]. A 
general discussion of its properties can be foimd in the books of Berger [9] and of GaUot, Hulin, 
Lafontaine [19]. 

Because of its importance, we resume here what we got in this special case. 

Proposition 4.14. Let pbea point of a smooth and connected Riemannian manifold (M, g), of dimension 

n, then the squared distance function from p is a locally semiconcave function on M and its gradient is an 
SBV vector field on M. Moreover, dp is C°° in M \ Cut{p) which is an open neighborhood of p. The cut 
locus ofp is C°° rectifiable with Hausdorff dimension at most (n — 1) and Cut{p) \ J is locally a finite 
union of smooth hypersurfaces. 

5. A Special Case 

In this section we study a special example which is interesting for the discussion of the next 
section. We assume that M is a two-dimensional analytic, connected and compact surface and K 
is an one-dimensional embedded analytic submanifold of M (a finite union of analytic curves). 
As before, our analysis also applies to closed sets K with analytic boundary. 
We look for topological results on the structure of the cut locus of K generalizing some arguments 
introduced to study the special case K = {p}, see Myers [25, 26]. Our goal is to show that 
Cut{K) is a finite graph and to connect its topological structure to the differential properties of 
the function . 

Clearly we have that UK, F and F^^, when it exists, are analytic. Notice that the fiber of UK is 
UpK^ {-1,1}. 

The strong result given by analyticity is that the number of the optimal focal points is finite. 

Lemma 5.1. The function c: UK R+ U {+00} defined in the previous section is analytic in the open 
set where it is not +00. 

Proof. With the same proof of Proposition 4.9, noticing that when n = 2 the set G coincides with 

Gi, we can show that the set of points (g, w, t) e UK x K.+ where dct dF(q, v,t) = is a finite 
union of analytic curves. Hence, as this set is the graph of the map c, we have the thesis. □ 

Proposition 5.2. The set J of optimal focal points of K is finite. 

Proof. Being M compact every geodesic carmot be minimal between its end points if it is longer 
than the diameter of M, hence a minimal geodesic joining an optimal focal point to K has to be 

shorter than a fixed constant. 

Consider an optimal focal point p and let F{q,v, t) be a minimal geodesic from K top which has a 
non invertible differential dF{q, v, a{q, v)) (notice that in this situation we have c{q, v) = cr{q, v)), 
we claim that ((j, v) is a critical point of the function c. 

By the Gauss Lemma (see [19, Chapter 3, Section E]), the differential dF{q, v, t) act on an element 
{w, s) G T(q^y^t)UK X M+ as follows, 

(5.1) dF{q,v,t){w,s) = duKF{q,v,t){w) + dtF{q,v,t){s) 

= X + sT 

where the two vectors X,T £ Tp(^q y i.^M are mutually orthogonal and T is the unit tangent 
vector to the geodesic F{q, v,t). Taking into accoimt that UK is locally a curve, this shows that if 
dF{q, V, t) is singular then duKF{q, v, t) = 0. 

Consider now the pull-back F*g of the metric tensor g on Tf^^ y i-^UK x IR+ via the map F. The set 
of points (g, v, t) where this form is not positive definite covers the graph of c. Computing this 
form using equation (5.1) we have, 

{F*9){q,v,t)i{w, s), {w, s)) = gF(q,v,t) {dF{q, V, t){w, s), dF{q, v, t){w, s)) 

(5-2) = + gF(q,v,t)iduKF{q,v,t){w), duKF{q,v,t){w)) 

= s'^ + h{q,v,t)gq{w,w) 
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for a non negative function h : UK x E+ R and where w G T(^q.v) UK is considered as a vector 
irxTgM. Clearly, the set of points where the function h : UK x R+ ^ R is equal to zero contains 
the graph of c by the previous discussion. 

Suppose that dc{q,v) ^ 0, recalling that UK is a curve there exists a small neighborhood B of 
{q, v) in UK where the map c is invertible and is analytic in the open set c{B) c M"*". Take a 
point (r, z) e B with c(r, z) < c{q, v) and consider the curve 7(5) : [0, c{q, v)] — > M defined by 

'7(0) = r 

< 7(s) = F{r, z, s) for s e (0, c(r, z)] , 
7(s) =F(c-^(s),s) for s e {c{r,z),c{q,v)] , 

that is, after a piece of geodesic, we follow the locus of first focal points. This is a piecewise ana- 
lytic curve in M starting from K and ending at the point p. The first piece of 7 is a geodesic, hence 
its length is c(r, z), the second piece follows the locus of first focal points. Using relation (5.2) we 
compute 

l^is)]^ = 9^is){dF{c-\s),s){w{s),l), dF{c-\s),s){w{s),l)) 
= 1 + h{c~'^{s),s)g{w{s),w{s)) 
= 1 

foreverys e {c{r, z),c{q,v)] and where «;(.s) = '^'■^^^ (.s). In the last equality we used the fact that 
the point (c~^ (s), s) belongs to the graph of c, where the function h is zero. 
Then the length of the second piece of 7 coincides with the variation in s, that is, c{q, v) - c(r, z). 
Finally the total length of 7 is c{q, v) = dxip)- 

Such curve is since the tangent vectors of its two parts are equal at the point F{r, z, c(r, z)), 
but it is not a geodesic for s G (c(r, z), c{q, v)], otherwise (by uniqueness) 7 should coincide with 
F{r, z, s) for every s € (0, c{q, v)] and this is impossible by construction. 

This fact implies that there must exist a shorter curve joining p with r and this is in contradiction 
with the assumption that F{q, v, t) is minimal, so the claim is proved. 

Arguing by contradiction, if the set of optimal focal points would be infinite then in a cormected 
component C of {c{q, v) < +00} c UK there would be infinite points (qi, Vi) e C where dc is 
zero and -F(gi, Vi,c{qi, Vi)) = pi are distinct optimal focal points. By compactness and the initial 
argument on the length of a minimal geodesic from an optimal focal point, there must exists an 
accumulation point of the set {{q^, Vi)} in C C {c{q, v) < L}, for a suitable constant L. Then, by 
the analyticity of c, this would imply that dc{q, v) is identically zero and c{q, v) is constant in the 
component C. 

Defining a function H : UK Ad by H{q, v) — F{q, v, c{q, v)) we have that 

dH{q, v) = duKF{q, v, c{q, v)) + dtF{q, v, c{q, v))dc{q, v) = , 
as duKF{q, v, c{q, v)) = and dc{q, v) = 0. 

So the map F{q, v, c{q, v)) is constant in C. This implies that all the points pi = F{qi, Vi, c{qi, Vi)) 
coincide contradicting the hj^otheses. □ 

As the map a : UK ^ M is continuous, the cut locus Cut{K) is given by a finite family of 
curves of kind s F{q{s), v{s), a{q{s),v{s))) where {q{s),v{s)) is a curve describing a connected 
component of UK. We say that p e Cut{K) is an end point if at the point p there arrives one and 
only one 1-cell of points of Cut{K). We are going to prove that every end point is an optimal 
focal point. 

First we exclude a very special case. 

Lemma 5.3. If at a point p e Cut{K) there arrive an infinite number of minimal geodesies then all these 
geodesies start from a unique connected component ofK which is a geodesic circle around p, that is the set 
of points of M at a certain distance Rfrom p. 

Moreover, p is an isolated point in Cut{K), more precisely Cut{K) n Bii{p) = {p}. Conversely, ifp is 
an isolated point in Cut{K) then there is a connected component of K which is a geodesic circle around p. 
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Remark 5.4. Notice that p is an optimal focal point by Lemma 4.8. 

Since the connected components of K are finite, it follows that the isolated point of Cut{K) are 
finite. 

Proof. If F{qi,Vi, t) is the infinite family of minimal geodesies F{qi,Vi, t) of length R ending at p, 
then all the distinct points qi belong to the geodesic circle of center p and radius Rin M. The 
set of points {{qi,Vi)} G UK clearly has an accumulation point, hence, by the analyticity of UK, 
the function H{q, v) = F{q, v, R) is constantly equal to p in the connected component of UK 
containing such accumulation point. Again by the analyticity of the connected components of 
UK and of the curves constituting K, we can conclude that the whole circle has to be a connected 
component of K. Hence, from every point of this circle there is a mirvimal geodesic ending at p 
and there cannot be other points of K inside the circle, otherwise their distance from p would be 
less than the radius R. 

Suppose now that p is isolated in Cut{K) and consider the open cormected component F of 

M \ K which contains p. The boundary of F is a subset K' of K and every minimal geodesic 
starting from K' with an initial velocity vector pointing toward F, must necessarily cease to be 
minimal at p, as Cut{K) n F = {p}. This last assertion follows from the fact that F, by point 6 of 
Proposition 4.6, can be continuously retracted on Cut{K) n F, hence Cut{K) n F is connected and 
then it coincides with {p}. This shows that there are irtfinite minimal geodesies from K top and 
we can conclude as in the first part of the lemma. □ 

Suppose now that p e Cut{K) \ J, so the number n > 1 of minimal geodesies F{qi, Vi,t) end- 
ing atp is finite. Consider a small ball B around p in M, then these n minimal geodesies cut the 
ball B inn sectors that we call S^. Any minimal geodesic starting in a sufficiently small neigh- 
borhood of {qi,Vi) e UK has its cut point in the ball B by continuity of the fimction a, moreover 
this geodesic cannot cross one of the geodesies F{qi,Vi,t) before reaching its cut point, otherwise 
this latter ceases to be minimal. Hence, considering the continuous curve of the cut points of 
the geodesies starting at the points of UK locally on the right side of {qi, Vi) (remember that U K 
is one-dimensional), we have that it is all contained in one of the sectors Si, more precisely, by 
continuity, in one of the two sectors adjacent to the geodesic F{qi, Vi,t). This curves gives a 1-cell 
of Cut{K) approaching p. 

With the same argument, considering the points locally on the left side of {qi,Vi) we obtain an- 
other 1-cell, in the other sector. 

Thus, we can conclude that the number of 1-cells of Cut{K) arriving at a point p is at least the 
number of the sectors Si, hence at least the number of minimal geodesies from pto K. 
This implies that every end point of Cut{K) where there arrives one and only one 1-cell, has a 
imique minimal geodesic to K so it has to be an optimal focal point. 

Putting together these facts and Lemma 5.3, by Proposition 5.2 the end points are finite. 
Following Myers [26], this result implies that the cut locus of is a linear graph and locally a 
tree, moreover the points where the order of the graph is greater than two are finite. 

Now we introduce the map from M to N U {oo} counting the number of minimal 

geodesies from K ioa point p. 

Proposition 5.5. An arc in Cut{K) containing no points of J and no interior points p with #Q{p) > 2 
is a regular analytic arc. 

Proof. Let 7 be such an arc in Cut{K). Consider a point pq E -y with F{qi,vi,(T{qi,vi)) = 
F{q2,V2,cr{q2,V2)) = Po, by the fact that po is not an optimal focal point, applying the im- 
plicit function theorem, there is an open neighborhood B of p in M \ K without optimal fo- 
cal points and there exist analytic functions zi,z2 : B UK, ti,t2 : B R+ such that, 
F{zi{p),tiip)) = F{z2{p),t2{p)) = p for every j9 G B and zi{B) n Z2{B) =0. If B is small 
enough, for every point p of B we have H=G{p) < 2, then ti{p) = t2{p) = dxip) if and only if 
p e Cut{K) n s = 7 n B. 

The rest of the proof proceed as in Proposition 4.7. □ 

Our last goal is to show that the order of a pointp G Cut{K), as a graph, is equal to #G{p). The 
order of a point p of Cut{K) is defined as the number of distinct 1-cells of Cut{K) arriving at p. 
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We have already seen before that that the order of p is always greater than the value #Q{p). 
We now prove the opposite inequality. 

Notice that the optimal geodesies cannot cross Cut{K), otherwise they cease to be minimal 
since they would intersect another minimal geodesic. We can take a small ball B aroxmd a point 
p e Cut{K) so that the 1-cells divide it in n sectors. Let S be one of these sectors and consider a 
sequence of points pi ^ Cut{K) all contained in S and converging to p such that F{qi ,Vi,t) are the 
minimal geodesies relative to Pi. By compactness, we can suppose that the points {qi,Vi) € UK 
converge to a point {q,v), hence the minimal geodesies F{qi,Vi,t) converge to a minimal geo- 
desic F{q, V, t) from K to p. Being the points pi contained in S the final part F{qi, Vi, t) H B oi 
the respective minimal geodesies have to be contained in the sector S and so also the final part 
F{q, v,t) (1 B oi the minimal geodesic for p. 

Taking into accoimt the fact that such minimal geodesic cannot intersect the cut locus, we con- 
clude that there is at least a minimal geodesic for every sector S. Being the number of the sectors 
equal to the order of Cut{K) as a graph, we proved the opposite inequality we claimed before. 
Hence, there are exactly =f^Q{p) 1-cells of the cut locus arriving at every point p e Cut{K). 
We summarize all the discussion of this section in the following theorem. 

Theorem 5.6. The set Cut{K) is a disjoint finite union of isolated points and linear graphs, each one 
locally a tree. 

The order of every point p G Cut{K) equals the function #G (p), counting the number of minimal geodesic 
from p to K. In particular, the set of points ofCut{K) with only one or more than two minimal geodesies 
is finite. 

The set of optimal focal points in Cut{K) is finite. 

All the isolated points and end points ofCut{K) are optimal focal points. 

Considering as vertices of the graph the optimal focal points and the points of order greater than two, the 
arcs connecting such vertices are regular analytic arcs. 

Remark 5.7. The analysis of this section also applies with small modifications to the case when 
K isa. finite set of points, considering an auxiliary set K consisting of a family of disjoint circles 
centered at the points of K with a radius R small enough. 

6. Singularities of Solutions of Hamilton-Jacobi Equations 

The ideas employed in the study of the distance from K and of the cut locus can be extended 
to analyse also the set of singularities of viscosity solutions of general Hamilton-Jacobi problems 



Moreover, geometric results on the cut locus suggest conjectures about the viscosity solutions of 
these equations. We give now an example. 

Suppose that A{x) is an analytic map from the closure of the imit ball B of to the space of 

positively defined 2 x 2-matrices. 
We consider the following problem. 



Using arguments similar to those of Section 5, it is possible to prove that the closure of the singu- 
lar set of the viscosity solution is a finite graph. 

If we look for the same result in the C°° case, a counterexample can be found as follows. It is pos- 
sible to endow the two-dimensional sphere with a C°° metric tensor g such that the cut locus of 
a certain point p is very wild, that is, it is not triangulable, hence it is not a finite graph (see [21]). 
Cutting away from the sphere a small geodesic disc D around p whose intersection with the 
cut locus of p is empty, and mapping stereographically from p the set § \ D on R^, we have that 
the closure of the singular set of the viscosity solution of Problem (6.2) in a ball of M^, where A 



(6.1) 




(6.2) 
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is given by the push-forward of the metric g via the stereographic projection, coincides with a 
homeomorphic image of the cut locus of p, hence it is not a finite graph. 

However, another results on the cut locus of a point says that for a generic C°° metric (in the 
category sense, the cut locus of every point of a surface is triangulable and has no points of order 
higher than three (see [30]). Hence, changing a little our point of view, this discussion suggests 
the following conjecture. 

Conjecture 6.1. For a generic C°° function A{x) from the closed unit ball B of to the space of 
positive definite 2 x 2-matrices, the closure of the singular set of the viscosity solution of problem 

{A{x)Vu{x), Vu{x)) = 1 inB , 
u = on dB 

is a finite graph. 

More in general, the same question can be asked about Problem 6.1 for a generic function H, 
domain fl or boimdary data uq. 
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